Quantum Simulation of Quantum Field Theories in Trapped Ions 
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We propose the quantum simulation of a fermion and an antifermion field modes interacting via 
a bosonic field mode, and present a possible implementation with two trapped ions. This quantum 
platform allows for the scalable add-up of bosonic and fermionic modes, and represents an avenue 
towards quantum simulations of quantum field theories in perturbative and nonperturbative regimes. 

PACS numbers: 03.67.Ac, 37.10.Ty, 03.70.+k 



A quantum simulator is a quantum device that can 
mimic the dynamics of another quantum system [I] [5] . 
It may represent a significant advance for studying the 
behavior of large quantum systems that cannot be sim- 
ulated efficiently with classical computers. Some pro- 
posals include the simulation of black holes in Bose- 
Einstein condensates |3], particle creation in expanding 
universes [H !Sj , quantum magnets [5] , and the emergence 
of relativistic physics and quantum field theories in ul- 
tracold atom lattices [Zl H]. In trapped-ion physics, a 
landmark was produced by the proposal for simulating 
the Dirac equation and Zitterhewegung in a single ion [5] , 
which was subsequently realized in the lab [T^. Recent 
advances in quantum simulations of relativistic quan- 
tum mechanics were the realization of the Klein para- 
dox [II1[T2], a proposal for implementing the Majorana 
equation and unphysical operations [13], and the quan- 
tum simulation of two Dirac particles interacting with a 
classical potential |14) . 

A natural step forward is to consider the constructive 
quantum simulation of quantized fields, in an effort to 
pave the way towards the implementation of quantum 
field theories (QFTs) in trapped-ion technologies. We 
point out that QFTs are significantly more involved and 
conceptually different from Refs. [9lfl4]. given that these 
previous works were related to relativistic quantum me- 
chanics, i.e., a single-particle theory. QFTs are among 
the most successful descriptions of the physical world, 
and give rise to the Standard Model of elementary parti- 
cles. There are different approaches to analize these the- 
ories, of which one of the most prominent is the Dyson 
series expansion in perturbation theory and Feynman di- 
agrams |15| . However, certain theories and parameter 
regimes cannot be analyzed in perturbation theory, like 
the strong coupling regime in which the coupling param- 
eter is so large that perturbative methods fail. Quantum 
chromodynamics, for example, cannot be studied in this 
way for low energies or long distances. New techniques 



have been developed for overcoming these difficulties, 
as lattice gauge theory computations |16| . A quantum 
simulator, on the other hand, could provide remarkable 
computational power to simulate QFTs faster than clas- 
sical computers. Among all plausible implementations, 
trapped ions may offer one of the most versatile and pow- 
erful, due to the high degree of quantum control [T7] . 

In this Letter, we propose a scalable quantum sim- 
ulation of interacting bosonic and fermionic quantum 
field modes in trapped ions. We consider first the case 
of a fermion and an antifermion field modes interact- 
ing via a bosonic field mode. This simplified model in- 
cludes already some of the interesting features of QFTs: 
particle creation and annihilation, self interactions and 
dressed states, and the possibility to study nonpertur- 
bative regimes in the ultrastrong coupling (USC) or 
deep strong coupling (DSC) regimes [18]. Moreover, we 
propose how to add progressively more fermionic and 
bosonic field modes with a scalable approach that may 
lead to the quantum simulation of a quantum field the- 
ory as quantum electrodynamics (QED). Actually, the 
use of lattice computations J16l to make nonperturba- 
tive assertions about a continuum quantum field theory 
requires a careful analysis. One has to consider critical 
points, correctly choosing lattice spacings and coupling 
constants. In this sense, it might well be that theoretical 
developments, akin to Wilson renormalization [19^, are 
needed so that our proposal can be used to make further 
predictions of quantum field theory phenomena. 

We will study a model under the following assump- 
tions: i) 1-1-1 dimensions; ii) scalar fermions and bosons; 
iii) one fermionic, one antifermionic, and one bosonic field 
modes. These considerations purport to make a realis- 
tic proposal with current trapped-ion experiments. Using 
conditions i) and ii) , we begin by considering a multimode 
bosonic and fermionic coupled system, with Hamiltonian 
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FIG. 1. Feynman diagrams for a number of terms in a finite- 
mode Dyson expansion associated with the interaction Hamil- 
tonian in Eq. (111. Note that loops do not imply sum over 
aU QFT modes. 
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(3) 
and Ofc are the bosonic annihilation operators that follow 
the commutation rules [afc,aj,,] = b(k — fc') while bp{dp) 
are fermionic(antifermionic) annihilation operators that 
obey the anticommutation rules {bp, b ,} = 6{p — p') and 

{dp,d ,} = 5{p — p'). Including spinors and polariza- 
tion amounts to incorporating additional modes and cou- 
plings, to be addressed with a many mode extension. 

In the interaction picture of Eqs. ^ and (Is]), the 
Hamiltonian reads, 



H = g dx ^\t,x)ip{t,x)A{t,x). 



(4) 



This is a simplified version of the QED Hamiltonian [T5] , 
allowing also the description of Yukawa interactions and 
the coupling of fermions to the Higgs field. 

Assumption iii) amounts to introducing two comov- 
ing fermionic and antifermionic modes. They will create 
single-excitation incoming wave packets when applied to 
the vacuum, at every time t in the Schrodinger picture. 
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Here, Qjj{pfj,p) are the fermionic (/) and antifermionic 
(/) envelopes of the incoming wave packets, that we 
consider as Gaussians centered in the average momenta, 
Pf f, with distant average initial positions. We intro- 
duce these propagating modes instead of the standard 
momentum eigenstates, given that the latter are delocal- 
ized over all space and are less realistic than the nor- 
malizable wave packets for describing physical particles. 
Although we deal with localized fermions, we consider 
a delocalized bosonic mode. In this way, we will have 
pair creation/annihilation processes and self interactions 
by considering just one bosonic field mode for the ease 
of experimental realization. We point out that our pro- 
posed quantum simulation of finite-number interacting 



quantized field modes includes all terms in a finite-mode 
Dyson expansion, instead of the standard one that takes 
into account all modes with a reduced number of per- 
turbative Feynman diagrams, see Fig. [11 Our motivation 
is to have a simplified model amenable to quantum sim- 
ulations with trapped ions, where the proposed theory 
emerges in the natural language of current experimental 
setups. The scalability of our approach may need specific 
developments of lattice theories [1^1 HO] . 

The fermionic modes b\^^ '^ ' id\^'^ ' , are the basis for de- 
scribing self-interacting dressed states, by emission and 
absorption of virtual bosons. They also represent the 
incoming states that will collide in a certain region of 
spacetime. These modes represent, at lowest order, the 
free evolution of the incoming wave packets. Here, the 
pair creation and annihilation is local, i.e., it takes place 
when the two wave packets of fermion and antifermion 
overlap, as corresponds to, e.g., the QED vertex. While, 
on the other hand, there is an interaction at a distance 
between fermion and antifermion mediated by the cou- 
pling to the bosons. 

Expanding the fermion field (pi in terms of these 
modes, and dropping the remaining anticommuting 
modes that complete the basis (not populated at lowest 
order), we get, in the interaction picture. 
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where the modes b[^ — j dp Qf{pf,p)bjj and dj^ 
J dp Qf{pjF,p)dj, do not depend on time. Here, 
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(9) 

are Fourier transforms, where we consider the ultrarela- 
tivistic limit, for simplicity, in which ojfj = \Pfj\- With 
this simplification, all time and space dependence in Eq. 
\4^ will be encoded in Qf{pf, x, t), Qjr{pjr, x, t), represent- 
ing Gaussian wave packets that propagate in spacetime. 



and in the time- and space-dependent phases that are 
associated to energy and momentum conservation. 
The bosonic field wih, in addition, be written as 

A{t, x) = aoe-*"°*e*'=°^ + ale"^°'e-''">'-' . (10) 

Consequently, the resulting interaction Hamiltonian is 

H = 9 Y. J''''iP^,P3^ko,t) epl^ ao + H.c, (11) 

where {O'^J.^fj- = &in,4i' and 
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(12) 



This Hamiltonian contains the self-interacting dynam- 
ics given by |/,/,n) O |/,/,n± 1) {\fj,n) denotes 
the state with one fermion, one antifermion, and n 



bosons), mediated by blJ)\-a_ak„ 
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It also includes pair creation and annihilation 



FIG. 2. (a) Setup for the trapped-ion simulation. (b) 
\{f,0,0\ip(t)i)\'^ as a function of t in units of too (red/upper 
curves), where \^p{t)i) is the evolved state from |i/;(0)i) = 
I/, 0,0), and average number of virtual bosons (blue/lower 

curves), (agao), for pi = O.lSwcO.lcJo, 0.05a;o,0.01a;o, 52 = 0. 
The largest amplitudes correspond to the largest couplings. 
(c) |(/, /, 0\il>{t)2)\^ as a function oit in units of ojq (red/upper 
left curve), where |'i/'(t)2) is the evolved state from 1-0(0)2) = 
I/, /,0), and average number of virtual bosons (blue/lower 
left curve), (ajao), for pi — 0.01a;o, 32 = 0.21tJo, (^t = 3/tJo, 
T = 30/ci;o, 5 = 0. (d) The same as (c) for gi = O.Icjq, 
32 = ^0, cTt = 4/ajo, T = 30/ioo, 5 = 0. 



processes given by \f,f,n) O |0, 0, n± 1), mediated by 



d\„bi 



in"in".s. 



and feJudinflfco i^. the quasi-resonant case, as well 



as dinbinOko 



and 6L4iaL 



map the 4-dimensional Hilbert space associated to the 
fermionic/antifermionic operators onto 4 internal levels 
of the first ion. For this, we consider a Jordan- Wigner 
mapping, bj^ = I ® a+ , bin = I ® (J~ , d\^ = cr+ ig) cr^ , 



and encode it in four internal levels of 

, |2), |3), |4), e.g., 6L = |4)(3| + |2)(1|, 



in the far off-resonant case. 
The last kind of transitions, as well as self-interactions, 
are off-resonant and would be neglected in the weak cou- 
pling regime, but would be allowed in our formalism for 
USC/DSC regimes [TH]- In our proposed setup, all per- 
turbative series terms are included, as shown in Fig. [Tl 

For practical purposes, we consider |fco| ^ wq, i-e., 
a slow massive boson. We may then approximate 
J^f-f{pf,Pf,ko,t) = J^f'f{pf,pf,ko,t) = giexp{-iujot), 
and J^f-f{pf,Pf,ko,t) ^ g2exp[-{t-T/2)^/{2<j^) + idt], ^ = -.9W ( |4)(1| aoc'"* + |1)(4| aje 



din = a~ (S) (J^, 
the first ion, |1) 



'^in = |4)(2| — |3)(1|, the vacuum state is state |1), and 
I/) = |2), I/) = -|3), I/,/) = -|4). Accordingly, Hamil- 
Ionian (131 results in 



where 52/51 gives the relative strength of the pair cre- 
ation with respect to the self- interaction, 6 = uif+ojf—uJo 
and T is the total time of the process, being at the tem- 
poral interval of the interaction region. Thus, the self- 
interactions are always on, while the pair creation and 
annihilation take place only when the fermion and an- 
tifermion wavepackets overlap, as they should. Accord- 
ingly, the Hamiltonian we aim to simulate is 



J „-M 



.g(i)( |1)(4| aoe-'(2c.o+5)t ^ |4^^;^| ^j^^^{2^o+S)t 
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(14) 



H = .9ie-'"°* [blhnao + dindLao) (13) 

e'^'bldlao + e-^t^-o+^^^din&inaoj + H.c. 



'526 



We propose to implement this Hamiltonian dynamics in 
a system of two trapped ions, see Fig. [2^. The bosonic 
mode will be encoded in the center-of-mass (COM) vi- 
bronic mode of the two-ion system. We envision to 



Here, the first line corresponds to a detuned red sideband 
interaction between |4) and |1) with time-dependent 
Rabi frequency g{t) = g2exp[-{t - T/2f/2a^]. The 
second line is a detuned blue sideband interaction, be- 
tween the same levels and with the same Rabi frequency. 
The third line can be developed applying detuned red 
and blue sideband interactions to |3) and |2) to pro- 
duce (|3)(2| — |2)(3|)[aoexp(— ia;ot)-l-aQexp(ia;oi)]/i, and 
a rotation of |3) and |2) with a classical field to pro- 



duce the change (|3)(2| - |2)(3|)/i — > (|3)(3| - |2)(2|). 
These operations are equivalent to a rotation from ay 
to az in the subspace of |3) and |2) and, therefore, 
not affecting the other Hamihonian terms. To imple- 
ment the last line we could just periodically drive the 
trap, but this method is only applicable to the COM 
mode. The alternative scalable approach that we pro- 
pose, applicable to any mode, consists of applying to 
the second ion detuned red and blue sidebands to gener- 
ate the term gi ax[aQexp{—iujQt) -I- aQexp(ia;oi)] where 
the matrix ax corresponds to the second ion internal 
levels. Preparing this ion into an eigenstate of a^, 
i.e (TxlV') = i|V')i we produce the desired interaction 
aoexp(— iwgt) -|- aJexp(iwoO [H]. Note that in the limit 
of large |i| in Eq. (14), we could also study the spectrum 



properties in a straightforward manner. 

We have numerically analyzed some of the features 
that could be simulated in the experiment: i) Fermion 
self-interaction-. In Fig. [2]d, we show the self-energy 
processes of a single fermion, which emits and reabsorbs 
bosons as it evolves. The average number of emitted 
bosons grows with the coupling. The effective Hamilto- 



nian that results from projecting Eq. ( 13 1 onto a single 

initial fermion, 

Ref. 



coincides with the DSC Hamiltonian of 
], Eq. (4), for qubit frequency equal to zero. Ac- 
cordingly, what is observed is emission and reabsorption 
of the bosons at a period which goes as 2tt/ujq, where 
Wo is in our case the simulated boson frequency, ii) 
Pair creation/annihilation-. In Fig. Pp, we show how 
a fermion/antifermion pair is annihilated, giving rise to 
a boson. Hi) Nonperturbative regime-. For large 52 val- 
ues, (72 > Wo, we enter into the nonperturbative regime, 
in which Feynman-diagram techniques are of little help. 
Here, we obtain that the number of created bosons is 
much larger due to the nonresonant terms (see Fig. [2a). 
The dynamics becomes more complex and strongly de- 
pendent on the specific coupling values. 

Having presented the basic three-mode model, we dis- 
cuss plausible steps in an effort to approach a full-fledged 
quantum field theory. We envision a feasible scalable 
method that allows a given trapped-ion setup to add-up 
more fermionic and bosonic field modes. The resulting 
Hamiltonian, equivalent to Eq. 
many modes, will be 



(11) with extension to 






^^32^1^ ^'''iP^,PJ,kl,t) 



'^l^l, ^k, 



H.C., (15) 



/ ij 



with i,j = l,...,n/, Uf the number of 
fermionic/antifermionic modes, and I = l,...,n;,, 
with rih the number of bosonic modes. Now the 
functions T'''^{pi,pj,ki,t) will contain the overlaps 
between the overcomplete set of Gaussians associated 
to the different incoming modes, Qi{pi,p), i = l,...,n/, 
similarly to Eqs. (12 1. We can consider additional field 
modes by involving more ions, including internal and 
external degrees of freedom, aiming at an intermediate 



scalability [H]. Here, it would be advantageous to make 
use of the transverse modes of a large ion string [22]. 
These modes can be closely spaced in frequency, such 
that it is not necessary to generate laser frequencies 
for each of them, relaxing infrastructural requirements 
while approximating the continuum to a certain extent. 
For the extension to many field modes, it is possible 
to consider the use of a digital simulator approach 
using many ions. The nonlocal spin operators could in 
principle be realized via recently developed stroboscopic 
techniques P5H?7] . We point out that already with 10 
two-level ions and 5 phonons per ion, one could perform 
quantum simulations of interacting quantum field modes 
that are beyond the reach of classical computers, that 
is, a Hilbert space dimension of 10^" ^ 2^^. 
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